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New universality class in percolation on multifractal scale-free planar stochastic lattice 
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We investigate site percolation on a weighted planar stochastic lattice (WPSL) which is a multi¬ 
fractal and whose dual is a scale-free network. Percolation is typically characterized by percolation 
threshold pc and by a set of critical exponents /3, 7 , n which describe the critical behavior of 
percolation probability P{p) ~ (pc — p)^, mean cluster size S ~ (pc — p)~^ and the correlation 
length ^ ~ (pc — p)”". Besides, the exponent r characterizes the cluster size distribution function 
ns{pc) ~ and the fractal dimension df the spanning cluster. We obtain an exact value for pc 
and for all these exponents. Our results suggest that the percolation on WPSL belong to a new 
universality class as its exponents do not share the same value as for all the existing planar lattices. 

PACS numbers: 61.43.Hv, 64.60.Ht, 68.03.Fg, 82.70.Dd 


Percolation is perhaps one of the most studied prob¬ 
lems in statistical physics because it provides a general 
framework of statistical theories that deal with structural 
and transport properties in porous or heterogeneous me¬ 
dia Bi- To study percolation one has to first choose a 
skeleton, namely an empty lattice or a graph. The model 
can then be defined by one sentence. Each site/bond of 
the lattice or graph is either occupied with probability p 
or remains empty with probability 1 — p independent of 
the state of its neighbors. For small values of p we see 
mostly single or a few contiguous occupied sites which 
are called clusters. As p increases, the mean cluster size 
always keeps growing at an increasingly faster rate till 
it comes to a state when suddenly a macroscopic cluster 
appears for the first time spanning from one end of the 
lattice to its opposite end. This sudden onset of a span¬ 
ning cluster in an infinite system occurs at a particular 
value of p known as the percolation threshold Pc- This is 
accompanied by sudden or abrupt change in the behav¬ 
ior of the observable quantities with a small change in its 
control parameter p. Such change is almost always found 
to be the signature of phase transition that occur in a 
wide range of phenomena ii- This is why scientists 
in general and physicists in particular find percolation 
theory so attractive. Indeed, the insight into the per¬ 
colation problem facilitates the understanding of phase 
transition and critical phenomena that lies at the heart 
of the modern development of statistical physics. 

Percolation on disordered lattice is potentially of great 
interest since many real-life phenomena deal with such 
disordered systems [^. In recent decades there has 
been a surge of research activities in studying percola¬ 
tion on random and scale-free network because the co¬ 
ordination number disorder of these networks is closely 
tied to many natural and man-made skeleton or medium 
through which percolation occurs. For instance, infec¬ 
tious diseases, computer viruses, pinion, rumors, etc 
spreads usually through networks [6l-[lll|. Besides, flow 
of fluids usually takes place through porous medium or 
through rocks and hence the architecture of the skele¬ 
ton is anything but regular ©[ii. In fact, transport 
of fluid through multifractal porous media such as sedi¬ 
mentary strata and in oil reservoir is of great interest in 


geological systems [l3, [T^ . In this rapid communication, 
we investigate percolation on weighted planar stochastic 
lattice (WPSL). One of us recently has shown that its co¬ 
ordination number distribution follows a power-law and 
its size distribution can be best described as multifractal 
[ 1 ^ . In contrast, scale-free networks too have power- 
law coordination number distribution but nodes or sites 
in the scale-free networks are neither embedded on spa¬ 
tial positions nor have edges or surfaces. Our goal is to 
find how the two aspects, multifractality and power-law 
coordination number distribution, leave their signature 
in the percolation processes. Classification of percola¬ 
tion into universality classes depending on the common 
sets of critical exponents has been of wide interest. It is 
well-known and widely accepted that the values of the 
critical exponents depend only on the dimension of the 
lattice and independent of its detailed structure and of 
the type of percolation, namely site or bond percolation. 
It is note worthy to mention that when a planar lattice 
is only multifractal but its dual is not a scale-free net¬ 
work, the resulting percolation still belongs to the same 
universality as the one for regular planar lattice El- We 
report for the first time that all the exponents for WPSL 
are completely different from the existing known values 
for the planar lattice {d = 2) and hence the percolation 
on WPSL belongs to a new universality class. 

The construction process of the WPSL starts with an 
initiator, say a square of unit area. The generator is then 
defined as the one that divides randomly first the ini¬ 
tiator into four smaller blocks. Thereafter at each step 
the generator is applied to only one block, by picking 
it preferentially with respect to areas, and divides them 
randomly into four smaller blocks. Each step of the di¬ 
vision process is defined as one time unit. It is thus a 
process of partitioning an initiator by random sequen¬ 
tial division that creates contiguous blocks of different 
sizes. The number of blocks N after time t therefore is 
N = 1 -k 3t and hence it grows albeit the sum of the areas 
of all the blocks remains the same. Thus, the number of 
blocks N increases with time at the expense of the size of 
the blocks. Indeed, the mean cell area decreases with N 
like (a) ^ N~^ or with t like (a) ^ since (a) = a/N 
where the total area a of all the blocks is always equal 
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to one. We therefore need to emphasize two things here. 
First, we can define the size of the side L of the WPSL 
as like we do for square lattice L = Second, 

to make the cells of the WPSL have the same size, in the 
statistical sense, as we increase N, we have to scale up 
the cell sizes by a factor of t. Else, we just need to multi¬ 
ply each quantity we measure by a factor t to compensate 
the decreasing size of the blocks. 

Perhaps, the construction process of WPSL is triv¬ 
ially simple but its various properties are far from simple. 
Firstly, it evolves following several non-trivial conserva¬ 
tion laws, namely ^ independent of time 

or size of the lattice V n, where Xi and yi are the length 
and width of the ith block. Secondly, its dual, obtained 
by replacing each block with a node at its center and com¬ 
mon border between blocks with an edge joining the two 
vertices, emerges as a scale-free network [^. Thirdly, if 
one considers that the fth block is populated with prob¬ 
ability Pi ^ xf or yf then the gth moment of pi can be 
shown to exhibit power-law Zq{S) ~ where S is the 

square root of the mean block area and 

T(g) = W + 16 - (3 q + 2). (1) 

Note that t(0) = 2 is the dimension of the WPSL and 
t(1) = 0 follows from the normalization of the proba¬ 
bilities X^Pi — 1 [ 3 . The Legendre transform of T(y) 
on the other hand gives the multifractal spectrum f(a) 
where the exponent a is the negative derive of T(q) with 
respect to q. Yet another features of the WPSL is that it 
emerges through evolution and the area size distribution 
function of its cells exhibits dynamic scaling [T^ . 

To study percolation on the WPSL we employ Ziff- 
Newman algorithm (2^ in which all the labeled sites or 
cells i = 1,2, 3,..., (1 -I- 3t) are first randomized and ar¬ 
ranged in an order in which the sites will be occupied. 
The good thing about this algorithm is that we can cre¬ 
ate percolation states consisting of n -I- 1 occupied sites 


simply by occupying one more site to its immediate past 
percolation state consisting of n occupied sites. Each 
time thereafter we occupy a site, it may happen that ei¬ 
ther an isolated cluster is formed or a group of contiguous 
sites linked by common border may get bigger either by 
agglomeration or by coagulation. We keep track of the 
number of clusters and their sizes as a function of n vis- 
a-vis the occupation probability p = n/{l -\- 3t). In fact, 
the product of the number of occupied sites n and the 
mean area 1/(1 -I- 3t) is equal to the mean area of all the 
occupied sites (a(n)) which is equal to 1 when all the 
sites are occupied i.e., when n = N. In our simulation 
we use periodic boundary condition where the lattice is 
viewed as a torus, thus without edge or surface, where 
sites are randomly occupied with probability p. 

In percolation, one of the primary objectives is to find 
the occupation probability Pc at which a cluster of con¬ 
tiguous occupied sites span the entire lattice, either hor¬ 
izontally or vertically, for the first time. Of course, the 
occupation probability at which it occurs at each inde¬ 
pendent realization on finite size lattice will not be the 
same. In reality, we can get spanning even at very much 
less than pc or not get it even at a much higher p than 
Pc- This is exactly why the percolation theory is a part 
of statistical physics. One way of dealing with this is to 
use the idea of spanning probability W(p) [2l| . Consider 
that we have performed m independent realizations and 
for each realization we check exactly at what value of 
p = n/N there appears a cluster that connects the two 
opposite ends either horizontally or vertically, whichever 
come first. The spanning probability W{p) is the proba¬ 
bility of occurrence of spanning cluster. It is obtained by 
finding the relative frequency of occurrence of spanning 
cluster out of m independent realizations. The plot in 
Fig. shows W{p) as a function of p for three dif¬ 
ferent lattice sizes. It thus represents the probability of 
finding a spanning cluster at occupation probability p for 
a fixed lattice size. One interesting point is that all the 
three plots meet at one particular point. It has a special 
significance as it means that if we could have data for in¬ 
finitely large lattice the resulting plot would also cross at 
the same meeting point. This meeting point is actually 
the percolation threshold Pc = 0.526846 for the WPSL. 

A careful look at the plots of Fig. ([5]) we find that if we 
increase L then a given fixed value of W is obtained at 
increasingly higher value of p for p < pc- To quantify this 
we draw a horizontal line, for instance at W{p) = 0.3, 
and a vertical line passing through the Pc value. Say, 
the horizontal line intersects all the three curves and the 
vertical line for different L at A, B, C and at O. We 
find that the distance OA, OB, OC etc which represents 
{pc — p) and plot them in the log-log scale as a function 
of t. The resulting plot gives a straight line with slope 
0.2966 ± 0.0055. Using L ~ we can write 

(Pc-p) ^ (2) 

where l/z/ ~ 0.6 01 v = 5/3. This is different and quite 
a bit higher than the known value z/ = 4/3 for all planar 
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Occupation probability p 


FIG. 2: Plots of spanning probability IF(p) vs occupation 
probability p for different lattice size. The vertical is drawn 
at Pc = 0.526846. In the inset, we plot W{p) vs (j> — Pc)L^^'' , 
where v = 5jS and find an excellent data collapse. 

lattices. For consistency check, one can now plot the p 
values at A,B,C etc versus . The intercept of the 
resulting linear fit gives the desired pc value and hence 
this offers an alternative method of measuring p^. The 
quantity ij)c—p)L^^'^ is a dimensionless quantity, accord¬ 
ing to Eq. in the sense that for a given value of W 
as L —>■ oo the value of (pc — p) —>■ 0 such that the nu¬ 
merical value of {pc—p)L^^’^ remains invariant regardless 
of the lattice size L. We now plot 'w{p) as a function of 
{pc-p)L^/'", see the inset of Fig. and find that all 
the distinct curves of Fig. ([3]) collapse onto a single uni¬ 
versal curve. It implies, according to finite size scaling 
hypothesis, that 

W{p)^L^cl>({p-p,)L^/''), (3) 

with exponent p = 0 where (j) is the scaling function . 
It states that the spanning probability W itself is a di¬ 
mensionless quantity provided it is measured in the scaled 
variable {pc — p)L^^'^ [1^. It also means that the span¬ 
ning probability for infinite lattice size would be like a 
step function around Pc- 

It is well-known that like Ising model percolation too 
display a continuous phase transition and hence like mag¬ 
netization of the Ising model there must be an order 
parameter for the percolation model too. The fact is 
that not all the occupied sites belong to spanning clus¬ 
ter. We thus can define the percolation probability P 
which must be zero below pc and should increase con¬ 
tinuously beyond Pc - & characteristic feature for order 
parameter. We define it as the ratio of the area of the 
spanning cluster Agpunt to the total area of the lattice at 
and hence P{p) = ^span since the the total area of the 
lattice is always equal to one. Unlike W{p) vs p the dis¬ 
tinct curves of the the P{p) vs p plots, see Fig. ([3]), for 
different size do not meet at one unique value, namely at 
Pc which we can only appreciate if we zoom in. Never¬ 
theless, following the same procedure we once again find 
(Pc —p)'^ with the same v value. Like for W{jp) if 

we plot P as a function (pc — p)P^^'' we do not get data 



FIG. 3: Plots of percolation probability P(p) vs p for three 
different size of the WPSL. In the inset we plot the same data 
but in the scaled variables and {pc — p)L^^'' and find 

an excellent data-collapse. 

collapse as before, instead we see that for a given value 
of {pc — p)L^^'^ the P value decreases with lattice size L 
following a power-law 

p ^ (4) 

where fijv = 0.135 ± 0.0076. It implies that for a given 
value of {pc the numerical value of PL^^" must 

remain invariant regardless of the lattice size of L. That 
is, if we now plot PL^/'' vs (pc — p)L^/'^ all the distinct 
plots of P vs p should collapse into a single universal 
curve. Indeed, such data-collapse is shown in the inset 
of Fig. ([3]) which implies that percolation probability P 
exhibits finite-size scaling 

P(pc - p, P) - L-°‘cj)(^{pc - . (5) 

Now, eliminating L in favor of pc — p in Eq. o we get 

P^{Vc-pf-, (6) 

where /3 ^ 0.225 or ,5 = 9/40 for WPSL whereas /3 = 
5/36 for all other known planar lattices. 

Percolation is all about clusters and hence the cluster 
size distribution function Ug (p) plays a central role in the 
description of the percolation theory. It is defined as the 
number of clusters of size s per site. The quantity sns(p) 
therefore is the probability that an arbitrary site belongs 
to a cluster of s sites and X]s=i i® probability that 
an arbitrary site belongs to a cluster of any size which is 
in fact equal to p. The mean cluster size S'(p) therefore 
is given by 

s(p)=y:»/. = £fP5i, (7) 

where the sum is over the finite clusters only. In the 
case of percolation on the WPSL, we regard s as the 
cluster area. It is important to mention that each time we 
evaluate the ratio of the second and the first moment of 
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(a) 


(b) 


FIG. 4: Mean cluster area S(p) as function of occupation 
probability p with different lattice size L. In the inset we plot 
the scaled variables SL vs (pc — and find that 

data for different system sizes are well collapsed in a single 
universal curve. 


of Ug we also have to multiply the result by t, the time at 
which the snapshot of the lattice is taken, to compensate 
the decreasing block size with increasing block number 
N. The mean cluster size therefore is 5” = i V x t 

where sug = p is the sum of the areas of all the 
clusters. Note that the spanning cluster is excluded from 
both the sums of Eq. ©• In Fig. 0 we plot S{p) as a 
function ofp for different lattice sizes L. We observe that 
there are two main effects as we increase the lattice size. 
First, we see that the mean cluster area always increases 
as we increase the occupation probability. However, as 
the p value approaches to Pc, we find that the peak height 
grows profoundly with L. 

The increase of the peak height can be quantified by 
plotting these heights as a function of L in the log-log 
scale and find 

5 ~ ( 8 ) 


FIG. 5: (a) The plot of log(ns) vs log(s) for different lattice 
size L. (b) The double-logarithmic plot of the size of the 
spanning cluster M against the lattice size L. 

where 7 = 2.883 which we can approximately write 
7 = 173/60. In contrast, 7 = 43/18 for all other pla¬ 
nar lattices. 

We can also obtain the exponent r by plotting the 
cluster area distribution function ng{pc) at Pc- We plot 
it in the log-log scale and find a straight line except near 
the tail. However, we also observe that as the lattice 
size increases the extent up to which we get a straight 
line having the same slope increases. It implies that if 
we performed on WPSL of infinitely large size we would 
have a perfect straight line obeying ng{pc) ~ s~'^ with 
T = 2.0724 which is less than its value for all other known 
planar lattices r = 187/91. We already know that mean 
cluster area S ^ 00 as p ^ Pc- According to Eq. ([7]), 
S can only diverge if its numerator diverges. Generally, 
we know that converges if a < — 1 and diverges 

if a > —1. Applying it into both numerator and denom¬ 
inator of Eq. o at Pc gives a bound that 2 < r < 3. 
Assuming 

ns(p) ~ ( 11 ) 


where "f jv = 1.73 ± 0.006321. A careful observation re¬ 
veals that there is also a shift in the p value at which the 
peaks occur. We find that the magnitude of this shift 
{Pc — p) becomes smaller with increasing L following a 
power-law {pc —p)'^ We now plot the same data 

in Fig (|4]) by measuring the mean cluster area S in unit 
of and {pc — p) in unit of respectively and find 

that all the distinct plots of S' vs p collapse into one uni¬ 
versal curve, see the inset of the same figure. It again 
implies that the mean cluster area too exhibits finite-size 
scaling 

S^lV((Pc-p)A'/"). (9) 

Eliminating L from Eq. (jS)) in favor of {pc — p) using 
(Pc — p) we find that the mean cluster area 

diverges 

( 10 ) 


and using it in Eq. © and taking continuum limit gives 

( 12 ) 

We know that sj diverges like (pc — pY^^ where a = 
Xjiydf') and hence comparing it with Eq. (jlOp we get 

r = 3 — 7 cr. (13) 

Note that the ramified nature of the spanning cluster 
at Pc is reminiscent of fractal. Indeed, we find that the 
the fractal dimension d/ of the spanning cluster can be 
obtained by finding the gradient of the plot of the size 
of the spanning cluster M as a function of lattice size 
L in the log-log scale (see Fig (IFbl) . We find dj = 1.865 
which can also be written as d/ = 373/200 for WPSL and 
that for square, triangular, honeycomb, Voronoi lattices 
is d/ = 93/48. Using the value of 7 and a in Eq. IT^ 
we get T = 2.0724 which we can approximately write as 
773/373. This is consistent with what we found from the 
slope of log[ns(pc)] vs log[s] plot shown in Fig. (15a)) . 


S^{pc-p)-^ 
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Exponents 

regular planar lattice 

WPSL 

u 

4/3 

5/3 

p 

5/36 

9/40 

7 

43/18 

173/60 

r 

187/91 

773/373 

df 

91/48 

373/200 


To summarize, we have studied percolation on a scale- 
free multifractal planar lattice. We obtained the Pc value 
and the characteristic exponents v, /3, 7 , r, a and df which 
characterize the percolation transition. Note that it is 
the sudden onset of a spanning cluster at the threshold 
Pc which is accompanied by discontinuity or divergence 
of some observable quantities at the threshold make the 
percolation transition a critical phenomena. One of the 
most interesting and useful aspects of percolation theory 
so far known is that the values of the various exponents 
depend only on the dimensionality of the lattice as they 
are found independent of the type of lattice (e.g., hexag¬ 
onal, triangular or square, etc.) and the type of percola¬ 
tion (site or bond). This central property of percolation 
theory is known as universality. Recently, Corso et al per¬ 


formed percolation on a particular mutifractal planar lat¬ 
tice whose coordination number distribution is, however, 
not scale-free like WPSL and still they found the expo¬ 
nents as for all the planar regular lattices [ 13 . Thus the 
most expected result would be to find a different value for 
Pc value as its coordination number distribution is totally 
different than any known planar lattice. However, finding 
a complete different set of values, see the table, for all the 
characteristic exponents was not expected since WPSL 
too a planar lattice. Interestingly, like existing values for 
regular planar lattices, the exponents of the values for 
WPSL too satisfy the scaling relations /3 = i'{d — df), 
7 = h'{2df — d), T = 1 + d/df. We can this conclude that 
percolation on WPSL belongs to a new universality class. 
It would be interesting to check the role of the exponents 
7 of the power-law coordination number distribution in 
the classification of universality classes. We intend to do 
it in our future endeavour. 
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reading of the manuscript and for valuable comments. 
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